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ABSTRACT 



Macieiewski et al.1 ()2010l ) claim to have detected a possible sinusoidal variation in 
the transit timing variations of exoplanet WASP-3b, which is currently the only known 
planet orbiting the star WASP-3. According to Maciejewski's analysis, this signal might 
be the consequence of gravitational perturbations caused by a hypothetical second ex- 
oplanet in the WASP-3 system. I report five transit timing measurements from the 
summer of 2010 which provide modest support for Maciejewski's proposed sinusoidal 
signal. 



1. Introduction 

Exoplanet transits are eclipses which occur when a planet orbiting another star passes directly 
between that star and the Earth. Although no current telescope can actually resolve the miniscule 
silhouette of the planet against its star, it is possible to detect the small fade produced as the planet 
blocks some of the light of its star from reaching Earth. Transits for a particular exoplanet occur 
at very regular intervals, and by examining a planet's transit timing variations (TTVs) — i.e., how 
early or late each transit occurs — it is theoretically possible to detect the effects of gravitational 
perturbations from other exoplanets in the same system. 



Macieiewski et al.l ()2010l ) observed six transits of exoplanet WASP-3b and combined them 
with eight others from previous publications. They detected a statistically significant sinusoidal 
signal in the fourteen TTVs, which spanned a period of over three years. Maciejewski found that 
the data are consistent with gravitational perturbations caused by a hypothetical second planet, 
WASP-3c, with about 15 Earth masses in a 2.02:1 orbital resonance with respect to WASP-3b - 
very close to a stable 2:1 resonance. They predict that such a planet would be difficult to detect 
via radial velocity measurements because the star's minute oscillations would obfuscate the weak 
radial velocity signature of WASP-3c. 

To test Maciejewski's claim, retired professional astronomer Bruce Gary calculated WASP-3b 
TTVs from dozens of observations in the Amateur Exoplanet Archive (AXA), his database of exo- 



- 2 - 



planet transit observations by amateur astronomers o He concluded that there was little evidence 
of sinusoidal variation in the amateur TTVs, many of which had relatively high uncertainties. 
Gary's procedure required each observer to perform basic data analysis, such as rejection of out- 
lier data and reference star selection, potentially injecting some measure of subjectivity into the 
process. Consequently, even though Gary performed the final analysis, it might not be appropriate 
to compare observations between different amateur observers, and his analysis is not necessarily 
conclusive. 



I observed five WASP-3b transits using an 11-inch Schmidt-Cassegrain telescope and a CCD 
camera at the observatory at the Jordan Hall of Science on the University of Notre Dame campus. 
Before each observing session, I synchronized the computer's time with Microsoft's Internet time 
server. I then took a continuous series of unfiltered images of WASP-3, spanning the length of the 
eclipse. When possible, I also obtained several hours of data before and after the transit in order 
to better characterize systematic errors. For example, the star's magnitude should theoretically 
remain constant outside of transit, but in practice, it will usually vary slightly as a function of the 
airmass if the target and reference stars are different colors. After each observing session, I reduced 
the data by applying flat fields and dark frames before performing aperture photometry to measure 
the fluxes of WASP-3 and twenty nearby reference stars. To account for night-to-night differences 
in the point-spread function of the stars, I used different aperture sizes for each session in an effort 
to maximize the quality of the photometry. 

In order to derive estimates of important transit parameters, especially the time of midtransit, 



WASP-3 AXA Light Curves, http://brucegary.net/AXA/WASP3/wasp3.htm1 



2. 



Observations 



Table 1. Transit Time Measurements 



Date 2010 



(UT) 



Mid- Transit 
(TBD) 




TTV 



Jun 4 
Jul 22 
Aug 2 
Aug 28 
Sep 8 



2455351.6832±0.0011 
2455399.6999±0.0015 
2455410.7802±0.0013 
2455436.6359±0.0008 
2455447.7155±0.0008 



0.2 ± 1.8 



0.9 ±2.2 



2.4 ± 1.6 



0.1 ± 1.2 



2.1 ± 1.1 
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I used an advanced spreadsheet designed by Gary for analyzing photometry of exoplanet tran- 
sits i Caltech's NASA Star and Exoplanet Database, which subsumed Gary's AXA in 2010, has 
specifically endorsed his procedure, which determines the best-fit model for each transit through 
X 2 -square minimization |£| This process produced estimates of a variety of important transit pa- 
rameters, including midtransit time, along with the corresponding uncertainties. To determine the 
expected time of midtransit, Tc, I adopted Maciejewski's ephemeris of 

T C {E) = 2454605. 56000 [BJD] + E x (1.8468355days) . 

The TTV is simply the difference between the predicted and observed times of midtransit. Table 1 
lists my individual TTV estimates using the above ephemeris. 

Since Maciejewski's ephemeris and subsequent TTV measurements relied upon the Barycentric 
Julian Date (BJD) reference frame and Terrestrial Time (TT) time standard, it was necessary to 
convert my data into a comparable time stamp. To accomplish this, I found the midtransit time 
in Julian Date and, with the assistance of an online conversion utility, I converted it to the BJD 
frame of reference and barycentric dynamical time standard (TBD)0 BJD_TT and BJD_TBD use 
the solar system's barycenter as their frame of reference, and according to Eastman, they vary 
from each other by no more th an 50 milliseconds, an insignificant difference for TTV analysis 



([Eastman. Siverd. &: Gaudill2010l ). These extremely accurate time stamps compensate for a variety 
of factors, such as Earth's orbital motion and gravitational perturbations from other planets in the 
solar system, which would otherwise taint TTV data due to the finite speed of light. 



3. Results 



I combined my data with Maciejewski's and fitted it with several different models using x 2 
minimization. A purely linear model was consistent with Maciejewski's period but the fit was 
quite poor, giving a x 2 parameter of 64.3 with 17 degrees of freedom. The resulting p-value of 
2.0 x 10~ 7 suggests an unacceptable fit to the data. I then fit a three-parameter sinusoid to the 
data, retaining the linear term because an inaccurate ephemeris would manifest itself as a linear 
variation in the TTV plot. I found that the combined data exhibited a sinusoidal variation with a 
period of 125.9±0.4 days with a full amplitude of 3.4±0.6 minutes. Furthermore, my analysis calls 
for a slightly refined ephemeris for WASP-3b, including a longer orbital period. Table 2 summarizes 
the parameters of my best-fit model, including my proposed ephemeris. 

Maciejewski proposes both short- and long-term sinusoidal variation in the TTV data, so I 



2 Bruce Gary's spreadsheets at |http: / /brucegary.net / book_EOA /xls.htm 

3 NASA/IPAC/NExScI, Star and Exoplanet Database, Summary Details for Amateur Light Curves, 



http://nsted.ipac.caltech.edu/NStED/docs/datasethelp/AXA.html 



'Jason Eastman's Barycentric Julian Date Converter, http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 
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Fig. 1. — The best-fit, single-periodicity model, with TTV estimates superimposed. The TTV is 
shown i n minutes. The yel l ow sy mbols show the new data presented here, while the other points 
are from lMacieiewski et al.l (|2010l ) and references therein. The revised ephermeris used in this plot 
is given in Table 2. The lower panel shows the residuals to the single-periodicity fit. 



performed the analysis again, allowing for two independent sinusoids. Without a linear term, the 
double sinusoid model produced a strong fit; however, once I included the requisite linear function 
to account for possible errors in the ephemeris, the reduced x 2 parameter plummeted from a nearly 
ideal 1.05 (x 2 = 13.62, 13 df) to just 0.5 (% 2 = 5.5, 11 df), indicating a serious overfit of the data. 
The extremely low reduced x 2 reflects that the total number of parameters in that model ballooned 
to eight, an ungainly number for just nineteen observations. Thus, I rejected this model, though I 
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was reluctant to do so because Maciejewski's computational models predicted that WASP-3c would 
cause two separate periodicities in WASP-3b's TTV. 



4. Discussion 

My observations provide equivocal support for Maciejewski's hypothesis. However, five closely- 
spaced observations with a small-aperture telescope are insufficient to provide convincing support 
for Maciejewski's sinusoidal model, which includes high-quality data obtained over a span of over 
three years. Additional observations would permit a more rigorous characterization of the TTV of 
WASP-3b, including an examination of whether there is a second periodicity in the TTV plot. 

One concern is that several of Maciejewski's TTV estimates between BJD 2454900 and 2455100 
are somewhat inconsistent with my fit. The residuals for these observations were much smaller for 
the double sinusoid model, but for reasons already discussed, I decided that a single periodicity 
best represented the data. Additionally, Maciejewski's published TTV estimate for BJD 2455102 is 
apparently incorrect. I used Maciejewski's ephemeris and estimated midtransit times to verify the 
published TTV estimates, and the BJD 2455102 datum showed the only discrepancy. For that date, 
Maciejewski's paper lists a TTV of 0.00120 days, but from other data in the paper I determined 
that the time of mid-transit corresponded to a TTV of 0.00058 days, a difference of 0.9 minutes. 
The most likely explanation is a misprint of that one TTV estimate in his table. 

Table 2. Best-Fit Model 



Full Amplitude (min) 3.7±0.6 

Sinusoid Period (days) 125.9±0.4 

Orbital Period (days) 1.8468372±0.0000006 

T (BJD) 2454605.55978±0.00010 

X 2 18.21 

df 14 

p-value 0.20 



The data quality for my 2010 June 4 observations left much to be desired, potentially making 
that TTV estimate unreliable. Jordan Hall is a notoriously vibration-prone building, and that 
night, high-frequency vibrations distorted the star images into elongated streaks. In some of the 
worst images, star images were over 13 arcseconds long. To put this into perspective, star images 
at Jordan Hall are normally circular and no larger than 5 arcseconds in diameter, even on a night 
of poor seeing. In light of this problem, it is prudent to be somewhat skeptical of the accuracy of 
this particular TTV measurement. 

On a final note, Maciejewski predicts that the hypothetical second planet would produce a 
transit depth of up to 0.35% if it undergoes transits. Although I did detect a barely significant, 
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0.3% fade in the flux of WASP-3 on 2010 August 1, my own follow-up observations strongly suggest 
that this feature was spurious, a product of a common systematic error, such as an imperfect flat- 
field. My experience suggests that while the equipment at Jordan Hall is just capable of detecting 
a transit depth of 0.3% under ideal conditions, it is extremely easy for any number of errors to 
produce false 'transits' of that depth. 

5. Conclusion 

Between June and September 2010, I observed five transits of WASP-3b, and my data, when 
analyzed in conjunction with previously published observations, suggests that there is indeed sig- 
nificant sinusoidal variation in WASP-3b's TTVs. Plainly, nineteen combined observations do not 
constitute a very large sample size, and more observations are necessary to explore this possibil- 
ity. Nevertheless, there appears to be considerable, albeit inconclusive, evidence that a second, 
undiscovered exoplanet is perturbing WASP-3b. 

I wish to thank Professor Peter Garnavich for his guidance with this project. Additionally, 
my research was subsidized in part by an AL/SCI UROP grant from Notre Dame's Institute for 
Scholarship in the Liberal Arts. 
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